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Abstract 

We discuss several new ideas for reactor neutrino oscillation experiments with a 
Large Liquid Scintillator Detector. We consider two different scenarios for a measure- 
ment of the small mixing angle #13 with a mobile v e source: a nuclear-powered ship, 
such as a submarine or an icebreaker, and a land-based scenario with a mobile reactor. 
The former setup can achieve a sensitivity to sin 2 2#i3 < 0.004 at the 90% confidence 
level, while the latter performs only slightly better than Double Chooz. Furthermore, 
we study the precision that can be achieved for the solar parameters, sin 2 2#i2 and 
Am|i , with a mobile reactor and with a conventional power station. With the mobile 
reactor, a precision slightly better than from current global fit data is possible, while 
with a power reactor, the accuracy can be reduced to less than 1%. Such a precision 
is crucial for testing theoretical models, e.g. quark-lepton complementarity. 
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1 Introduction 



Reactor experiments have always played a crucial role in neutrino physics. The first experi- 
mental evidence for the existence of the neutrino came from a reactor experiment [1,2] and 
the most precise measurement of the oscillation parameter Am^ has been performed by the 
KamLAND experiment in Japan [3], where the z/ e -disappearance of neutrinos coming from 
the surrounding power plants at an average distance of approximately 180 km has been 
studied. KamLAND data has uniquely identified the large mixing angle scenario as the 
correct solution of the solar neutrino problem. Furthermore, the non-observation of reactor 
neutrino disappearance at baselines ~ 1 km in the CHOOZ experiment in France is crucial 
for the current upper bound on the small mixing angle sin 2 2#i3 < 0.1 [4]. The role of future 
reactor neutrino experiments in measuring or constraining the value of #13 has already been 
extensively studied in the literature [5-8]. Usually it is proposed to construct a relatively 
small detector close to a large nuclear power station. However, it is important to keep in 
mind that the performance of such a measurement besides the systematical uncertainties 
depends only on the total exposure, which is proportional to the product of the thermal 
power of the reactor, the detector mass, and the running time of the experiment. Thus, it 
is also imaginable to take advantage of a bigger multi-purpose detector, such as for example 
the Large Liquid Scintillator Detector LENA [9,10], and combine it with a relatively small 
reactor. In this work, we consider using a mobile reactor, which will lead to interesting can- 
cellations of systematical errors. Mobile reactors are widely used on nuclear-powered ships 
or submarines, but there are also ideas to construct land-based removable power stations 
that may be used to deliver electricity to very remote areas [11]. 

Apart from searching for nonzero #13, an advanced reactor neutrino experiment can also 
provide a precision measurement of the solar oscillation parameters. It could improve the 
bounds on Am^, which are currently dominated by KamLAND, as well as those on the 
oscillation amplitude sin 2 2# 12 , which are currently dominated by solar data. This improved 
precision is required, for example, to test quark-lepton complementarity [12] to a high preci- 
sion, and to determine the neutrino mass hierarchy in 0z//3/3-experiments [13,14] and through 
the day-night effect in solar neutrino data [15]. A possible setup for such an advanced re- 
actor neutrino experiment aiming at a precision measurement of the solar parameters has 
already been presented in [16]. We will complement the discussion given there by studying a 
setup involving a rather small, possibly mobile, nuclear reactor, and a LENA-like detector. 
Furthermore we will discuss the possibility of constructing such a large detector close to a 
nuclear power station. Since it is usually not desirable to construct a Large Liquid Scintil- 
lator Detector close to such a power plant since it would then be blind to Geo-neutrinos, 
supernova relic neutrinos, and other very faint sources, one would have to choose a reactor 
which is scheduled for permanent shutdown after several years of data taking, or a newly 
built power plant which only starts operation after the low-background measurements have 
been completed. 

The outline of the paper is as follows: In Section 2, we will briefly review the neutrino 
oscillation framework and the underlying phenomenology of z> e -disappearance measurements 
in reactor experiments. In Section 3, we will discuss the prospects of 9 13 measurements 
with mobile reactors for two different scenarios, a nuclear-powered ship and a land-based 
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mobile nuclear reactor. We will first describe the expectations from analytical estimates, and 
thereafter present results based on numerical calculations. We carefully include systematical 
uncertainties, which are the main limitation to the achievable sensitivity. Next, in Section 4, 
we will discuss precision measurements of the solar parameters sin 2 20 i2 and Am^ with both, 
a mobile reactor and a large nuclear power plant. Again, we will give analytical estimates 
as well as numerical results. The latter will include a careful treatment of the Geo-neutrino 
background which turns out to have a significant impact on the optimization of the baseline. 

2 The neutrino oscillation framework 

The relevant oscillation channel for a reactor neutrino experiment is v e — > v e disappearance. 
The corresponding exact formula for the vacuum survival probability of an electron anti- 
neutrino V e with energy E at a baseline L is given by [17]: 

P{p e -> P c ) = 1 - sin 2 26>i3 sin 2 A 31 - (cos 4 13 sin 2 20 12 + sin 2 12 sin 2 20 13 ) sin 2 A 21 + 

I 

+ sin 2 6>i2 sin 2 20 i3 (- sin 2A 2i sin 2A 31 + 2 sin 2 A 31 sin 2 A 2 i), (1) 

where A 2 i = Am 21 L/ (4i?) and A 3i = Am 31 L/ (4_E). The parameters involved are the "reac- 
tor angle" 0i 3 , the "atmospheric mass squared difference" Am 31 , and the "solar parameters" 
6>i 2 and Am 21 . The solar parameters have been measured by the SNO and KamLAND ex- 
periments [3,18-21], and due to the observed MSW effect [22,23], the sign of Am 21 is also 
known. The atmospheric parameters have been measured by Super-Kamiokande [24,25], 
K2K [26,27], and MINOS [28,29], but the sign of Am 2 x has not been determined yet. Con- 
sidering the small mixing angle 0i 3 , there exists only an upper bound which is dominated 
by the CHOOZ experiment [4]. The best-fit values and allowed intervals for all oscillation 
parameters from a global three-flavor analysis can be found in [30-33]. 

While our numerical simulations are based on the full three-flavor treatment, we will only 
use zeroth-order approximations for the analytical estimates. These approximations are 
given by 

PsJpe -> P e ) » 1 - sin 2 20i 2 sin 2 A 2 l, (2) 
for the measurement of the solar parameters, and by 

^PatmOe -> v e ) « 1 - sin 2 20 i3 sin 2 A 31 (3) 

for the (13)-oscillations and are justified by the numerical values of the parameters. Since 
|Am 31 | ^> Am 21 , the (13)-oscillation length is much smaller than the (12)-oscillation length, 
so a typical 0i 3 experiment located at the first atmospheric oscillation maximum is hardly 
affected by the solar terms. On the other hand, for the amplitudes it holds that sin 2 20 i2 ^> 
sin 2 20i 3 , so the (13)-oscillation constitutes a negligible perturbation in a measurement of 
012 and Am 21 . 
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For the oscillation parameters, we assume the following true values in the numerical simu- 
lations: 

Am 2 31 = 2.5-10" 3 eV 2 , sin 2 2# 23 = 1, () 

Amli = 8.2-10" 5 eV 2 , sin 2 26 12 = 0.83. W 

3 0i3 measurement with a mobile reactor 

In this section we consider the potential of a mobile electron anti-neutrino source, combined 
with a Large Liquid Scintillator Detector, to measure the small mixing angle #13. We 
assume that the reactor is placed at two different baselines, the "near" and "far" positions, 
consecutively. The major limitation for such an experiment will be the systematical errors. 
If these are not under control, any deficit in the observed neutrino flux could be attributed 
either to neutrino oscillations or to a systematical bias in the initial reactor neutrino flux 
or in the detector properties. We will consider two scenarios with different systematical 
uncertainties: 

• A nuclear-powered ship, e.g. a submarine or an icebreaker (c/. [34]): 

Since the reactor is not required to be shut down in order to change its position, and 
since it can be moved between the near and far positions frequently, it is reasonable 
to assume that the unoscillated neutrino flux and spectrum are the same at both 
positions. This will turn out to be crucial for the cancellation of the associated sys- 
tematical errors. Such a nuclear-ship scenario could only be realized, if the detector 
were located under water, e.g. in the Mediterranean Sea. 

• A land-based nuclear reactor, such as the SSTAR design [11]: 

Here, the reactor needs to be shut down in order to be maneuverable, and such a 
movement will only be done once, so the neutrino flux and spectrum at the near and 
far positions will be uncorrelated. This problem could be alleviated by using a small 
dedicated near detector to reduce the uncertainties to those associated with the near 
detector, which are typically smaller than those associated with the reactor. 

For both scenarios, we consider a detector with a fiducial mass of 45 kt, which corresponds 
to the proposed LENA detector [9,10]. Note that some of our conclusions can also be applied 
to reactor neutrino experiments in which not the source, but the detector is mobile [35]. 
In such a scenario, the detector will only be moved once, so it is similar to our land-based 
scenario. However, the detailed impact of the systematical errors is different: If the reactor 
is movable, only the neutrino flux and spectrum uncertainties will be uncorrelated in the 
two phases of the experiment, while for a movable detector scenario, also at least some of the 
detector-related uncertainties should remain uncorrelated. Apart from this, the feasibility is 
different since mobile nuclear reactors do already exist and their development is continuing, 
while there is no practical experience on mobile neutrino detectors yet. Note, that only for 
the nuclear-powered ship scenario it can be assumed that all neutrino flux and spectrum 
associated systematical uncertainties are completely correlated in the two phases of the 
experiment, since, as mentioned above, a frequent change between the near and far position 
is possible. 
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To study the sensitivity of reactor neutrino experiments, we use a x 2 analysis, incorporating 
pull terms for the proper treatment of systematical uncertainties. In our numerical calcula- 
tions we assume the events to follow a Poisson distribution, but for illustrative purposes it is 
sufficient to consider the Gaussian approximation, which is very good due to the large event 
rates in a #13 reactor experiment. For the nuclear-powered ship scenario, our \ 2 expression 
has the form 



Here, Nf and Nf denote the event numbers in the i-th bin at the near and far positions, 
respectively. These event rates are calculated with the values for the oscillation parameters 
given in Eq. (j3J). Correspondingly, Tf are the theoretically predicted event rates for certain 
fit values of the oscillation parameters, and for the systematical biases aj and b. The second 
line of Eq. (jSJ) contains the pull terms which represent prior knowledge about the system- 
atics parameters. They give a penalty to biases that are much larger than the estimated 
systematical errors. In detail, we introduce the following systematical error sources: 

• The flux normalization uncertainty cr nor m = 2.0%, which is correlated between the 
near and far positions in the nuclear-powered ship scenario, but not in the land-based 
scenario. It accounts for the limited accuracy with which the thermal reactor power 
and thus the emitted neutrino flux can be determined. 

• The detector normalization error cidet = 2.0%, which describes uncertainties associated 
with the fiducial detector mass, the cross sections, the scintillator properties, and the 
analysis cuts. 

• A background flux error o"bk g = 10.0%. As backgrounds, we take into account Geo- 
neutrinos, amounting to about 1450 events per year [36-38], and a diffuse reactor 
background of about 850 events per year, which corresponds to the estimated flux 
from the 20 closest reactors at the proposed LENA site in Pyhasalmi (Finland) [39,40]. 
By taking into account all nuclear reactors in the world, the background rate would 
only increase marginally. 

• The shape uncertainty, o" s h ap e = 2.0%, which describes the limited knowledge of the re- 
actor neutrino spectrum. The corresponding parameter a s h ap e,i is independent for each 
bin i, hence the index i is introduced. This parameterization follows the discussion 



• An energy calibration error = 2.0%, which is parameterized by b. 

These errors are summarized in Table [H The dependence of the Tf on the systematical 
biases is given by 




A=N,F i 





in [6]. 
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Reactor Neutrino Flux 2.0% 

Detector Normalization 2.0% 

Detector Energy Calibration 2.0% 

Shape Error 2.0% per bin 

Normalization of 1% background 10.0% 

Table 1: The systematical uncertainties that are assumed in our numerical simulations (based on [8,41]). 

Tf- = (1 + a norm + a dct + a ehapeji )Sf(b) + (1 + a bkg + a dct )Bf(b), (6) 

where in turn Sf(b) nd Bf(b) are the signal and background rates for the wrong energy 
binning implied by nonzero b. They are obtained from the correctly binned rates Sf and 
Bf according to 

Sf(b) = (1 + b)[(S [si+1 - S ls} )(5 - [5\) + S [6i ], (7) 
5 = b-(i + t + (8) 

and a similar expression for Bf (see also ref. [6]). The quantity to m Eq. ([8]) is the energy 
threshold of the detector, expressed in units of the bin width. We have used the Gauss 
bracket |_ - J to denote the floor function. The expression in square brackets in Eq. ([7]) is 
essentially a linear interpolation between the events in bin |_<5J and those in bin \_5\ + 1. If b 
is not too large, the energy calibration never changes by more than the bin width, so that 
[5\ = i. The factor (1 + 6) in front accounts for the change of the bin width implied by b. 

The x 2 expression for the land-based scenario is similar to Eq. §5§ , but since the reactor flux 
and spectrum are uncorrelated between the near and far positions in this scenario, a nor m 
and a s ha P e,i get an additional index A = N, F. 

3.1 Analytical estimates 

We will now show that, for the nuclear-powered ship scenario, the most important system- 
atical errors can be eliminated in the actual measurement of #13 if the reactor is positioned 
at two different baselines consecutively. To quantify the performance of the experiment, 
we consider the sensitivity to sin 2 2^13, which is defined as the limit that can be set on 
sin 2 2#i3 assuming that the true value is 0. This quantity can be calculated by comparing 
the simulated event rates for non-zero test values of sin 2 2#i3 with those for sin 2 26*13 = in 
a x 2 analysis. 

For an analytical estimate, we consider a simplified version of Eq. ([3]), keeping only the 
reactor flux error <7 norm and the detector normalization error o"d e t, but neglecting spectral 
uncertainties, energy calibration errors, and backgrounds: 



* 2 = E E^[^(l + a norm + a dc 0(l-sin 2 2^ 13 sin 2 Af)-iV^%%^ + ^. (9) 



A=N,F i=l 



norm det 



Here, the index i runs over all n energy bins, while A takes the values N and F for the 
near and far baselines, respectively. Nf is the total event rate in the i-th bin at baseline L A 
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without oscillations, and Af = Am| 1 L j4 /4£'j is the oscillation phase for baseline L A and 
energy 

Since the systematical biases a nor m and eidet as well as the oscillation amplitude sin 2 26*13 are 
small, Eq. ^ can be approximated by 

n 22 
* 2 = E E^K°™ + a ^-sin 2 20 13 sin 2 Af) 2 + %2L + ^. (10) 

A = ]\[ p j = i ^norm °det 

flnorm and a Act are fitted by minimizing x 2 , therefore we calculate 

E E 2iV ^ + a ^ - sin2 2 ^3 sm 2 Af ) + ^= (11) 



C&norm A=7V _F 2=1 ^norm 



TT^ = E E 2N ? + - sin 2 29 13 sin 2 Af ) + . (12) 

By requiring the expressions fill I) and (j!2p to be zero and taking their difference, we obtain 



2 

^norm 
CT dct 



Cldet- (13) 



Substituting this back into Eq. (TTTj) or (Tl2|) , and assuming the iV^ 4 to be very large, the 
systematical biases can be expressed as 

sin 2 20 13 E,wEILi^sin 2 Af 
adet 1 4. n 2 ifji v V™ AP 1 ' ^ ' 

1 ' "norm/ u det Z^A=iV ,F Z^i=l iv j 



sin 2 2^3 EA=^Er=i^sin 2 Af 



-norm ~ i , 2 / 2 V- V" /V* A 

1 ^ °det/ "norm 2^A=N,F Z^i=l ly i 

Thus, in the limit Nf — > oo, x 2 is approximately given by 



(15) 



/V V™ 7V B sin 2 A B \ 

2 • 4 /) \^. r 4 Z^B=N,F 2^j=l iv j b m X 

X =sm 2^i3 2^ N i \ — ajB sm A i ■ ( 16 ) 

A=N,F i=l \ 2^B=N,F 2^3=1^3 J 



The la sensitivity limit for sin 2# 13 is determined by the condition \ — 1 ; 

n 



(jfsin 2 26 



13/ 



i v 

A=N,F i=l 



/ Eb=n,f EU Nf [sin 2 Af - sin 2 Af ^ ' 



iV/^oo 

i— > 0. 



l_^B=N,F Z^j=l iv j 

(17) 

This means that the sensitivity will be good if Nf is large and | sin 2 Af — sin 2 Af\ is not 
too close to zero. Therefore an optimal experiment involves one measurement very close to 
the reactor (sin Af » 0) and one around the first oscillation maximum (sin 2 Af « 1). 

Now compare this to the case where the reactor position is fixed, i.e. L N = L F . Then, 
= N[ = Ni/2 and sin 2 A^ = sin 2 Af = sin 2 A^. If we perform only a total rate analysis 
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(n = 1), the right hand side of Eq. (fl6l) will be very close to zero. This means that we have 
to consider the next-order term: 



2 _ Q norm u det _ -4 n I Z^i 

X - -5 h - sin 2#i3 ^ 77 — : — 2 



1IAV1 I 1 

/ V, iVi sin 2 A 



norm 



In this case, cr(sin 2 2^i 3 ) is of order a/ o"det + °norm- 

In reality, n will be chosen greater than 1 in order to exploit spectral information. In this 
case, the right hand side of Eq. ( fl6l) can be interpreted as a comparison of the oscillation 
phase in bin i, sin 2 Af, to the average oscillation phase Eb Ej Nf sin 2 / J2b Ej Nf . 
Therefore, this term is sensitive to the spectral distortion caused by neutrino oscillations, 
and in principle, the systematical errors cr norm and o"d e t can be eliminated without going to 
different baselines. This seems reasonable because the L/E dependence of the oscillation 
probability Eq. ([3]) implies that measuring at different energies is equivalent to measuring 
at different baselines. However, our numerical results presented below will show, that in 
reality it is still highly advantageous to change the baseline. 



3.2 Numerical results 



To demonstrate the full impact of systematical errors on the sin 2 2#i3 sensitivity, we have 
performed numerical simulations with a modified version of the GLoBES software [42], which 
incorporates a \ 2 expression similar to Eq. (jSJ), but is based on the more realistic assumption 
of Poisson statistics and uses the exact three-flavor oscillation probability without applying 
any approximations. The reactor spectrum has been taken from [3,43], and the cross 
sections are from [44]. We assume an energy resolution of 0.091 ^/E/MeV MeV ( [10,45]) 
and sort the events into 67 energy bins. Additionally, we assume that the reactor is located 
at baseline L\ for i tot L 2 / (L\ + L 2 ) and at baseline L 2 for t tot L 2 / [L\ + L|), where t tot is the 
total running time of the experiment. This ensures that the numbers of events in the near 
and far detectors are comparable. We take L\ = 0.1 km and L 2 = 1.3 km. 

Figure [1] shows the sensitivity of different scenarios as a function of the total exposure 
$. In addition to the nuclear-powered ship and land-based scenarios that where discussed 
above, we also show for comparison the performance of a setup similar to Double Chooz [8] 
with one reactor and two detectors. In this Double Chooz-like scenario, the uncorrelated 
detector normalization and energy calibration errors are taken to be only 0.6% and 0.5% 
because many detector-side effects will cancel. For very low exposures, $ < 0.02 GW kt y, 
all three scenarios are limited by the overall statistics and by backgrounds. As long as 
bin-to-bin errors, that are uncorrelated between the different bins as a s h a pe and additionally 
uncorrelated between the near and far phase, are neglected (see left plot in Fig. [[]), the 
nuclear-powered ship scenario (solid blue/black curve) follows the statistical limit also for 
larger exposures because systematical errors cancel completely in this scenario, as expected 
from our analytical estimates. For the land-based scenario (solid cyan/gray curve), only the 
systematical errors associated with the detector cancel because we have assumed the reactor 
neutrino flux to be completely uncorrelated before and after the displacement of the reactor. 
For the spectrum, we assume a partial correlation, which is parameterized by introducing in 
addition to the shape error from Eq. (jSJ) a new bin-dependent bias with a 1% error, which is 
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Figure 1: Sensitivity to sin 2 29i3 at the 90% confidence level as a function of the exposure for a nuclear- 
powered ship scenario (solid blue/black curve), a land-based mobile reactor scenario (solid cyan/grey curve), 
a land-based mobile reactor scenario with an associated near detector (dashed cyan/grey curve), and a 
Double Chooz-like setup (dotted red/black curve) with one reactor and near and far detectors. The left 
plot has been calculated under the assumption of no bin-to-bin errors, while the right one includes an 
uncertainty (Thin-bin = 0.5%. 

only present at baseline L2. As soon as the exposure exceeds 0.02 GW kt y, the sensitivity 
of the land-based scenario becomes limited by the uncorrelated flux normalization error. As 
the exposure increases further, this error becomes less dominant because the event numbers 
are so large that spectral information can be exploited. However, above 10 GW kt y, the 
uncorrelated shape error prevents a further increase in the sensitivity. This can be avoided 
if a small near detector with a mass of 45 t, an uncorrelated normalization error of 0.6%, 
and an uncorrelated energy calibration error of 0.5% is employed (dashed cyan/gray curve). 
Finally, in the Double Chooz-like scenario (dotted red/black curve), the errors associated 
with the reactor cancel, but those associated with the detectors remain. As the reactor 
flux has larger uncertainties than the detector normalization, this scenario is better than 
the land-based mobile reactor scenario for luminosities of 0.02 GW kt y - 6 GW kt y. At 
the onset of the spectrally dominated regime, the two curves meet again, but for very large 
luminosities the Double Chooz-like near/far scenario is again better than the land-based 
scenario (without associated near detector) because the shape error is canceled by the near 
detector. 

The vertical lines in Figure [T] indicate the total exposures that are to be expected from 
the Double Chooz experiment, its possible upgrade scenario Triple Chooz (see [41]), and 
a possible total exposure of the mobile scenarios (0.1 GW thermal power of the mobile 
reactor source, 3 years data taking, 45 kt fiducial mass of the LENA detector). While the 
performance of the land-based mobile reactor scenario is comparable to the Double Chooz- 
like scenario (dotted red/black curve) taken at the same total exposure, a nuclear-powered 
ship scenario could, due to the excellent cancellation of systematical uncertainties, achieve 
a sensitivity limit sin 2 2#i3 < 0.003 at the 90% confidence level. 
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Sensitivity to sin 2 26» I3 at 90% C.L. Sensitivity to sin 2 20 13 at 90% C.L. 
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Figure 2: Sensitivity limit to sin 2 2#i3 at the 90% confidence level for the nuclear-powered ship scenario 
(left) and the land-based scenario (right) as function of the near and far baselines. 

Even under the assumption of an uncorrelated bin-to-bin error o~bm-bin = 0.5% (see right 
hand plot of Fig [I]), this excellent sensitivity decreases only marginally. crbin-bin parame- 
terizes all kinds of unknown backgrounds and detector non-linearities. It is introduced in 
Eq. (jSJ) in a similar way as cr shape , but the corresponding parameter a£ in _ hini depends not 
only on the bin, but also on the detector position, hence the index A = N, F for near and 
far phase is introduced. 

Let us now address the issue of choosing the optimal baselines for the two phases of the 
experiment. Figure [2] shows that, for the nuclear-powered ship scenario, it is clearly advan- 
tageous to choose L\ very close to the reactor in order to measure the unoscillated flux (0th 
oscillation minimum), and L2 around the first oscillation maximum at about 1.3 km. For 
the land-based scenario, the overall sensitivity is worse, but now the baseline combination 
1st minimum/lst maximum gives better results than the combination 0th minimum/lst 
maximum. From a practical point of view, this turns out to be advantageous because in the 
land-based scenario it will usually not be possible to move the reactor close to the detector, 
as the latter will be located deep underground. We have also studied scenarios in which 
more than two different baselines are used. However, no significant improvement in the 
sensitivity could be achieved this way. 

4 Measurement of the solar parameters 

An often neglected topic is the issue of a precise measurement of the so-called solar oscillation 
parameters (9 12 and Am^) with a reactor neutrino experiment located at the first solar 
oscillation maximum (earlier studies can be found in [16,46,47]). Although e.g. the value 
of #12 has been determined by solar experiments and KamLAND [33] to be about 33° ± 3°, 
which seems to be quite accurate, it should be stressed that the relative uncertainty is still 
about 10%. Since precise knowledge of 9 12 is very important, in particular for the distinction 
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between the normal and inverted neutrino mass hierarchy in Oz//3/3-experiments [13,14] and 
via the day-night effect of solar neutrinos [15], as well as for exploring the possible existence 
of quark-lepton complementarity [12], a more precise measurement is desirable, and we will 
show that it can be performed with a Large Liquid Scintillator Detector. 

For the measurement of the solar oscillation parameters, two different scenarios are consid- 
ered: the first one (" SMALL") employs a mobile reactor with a thermal power of 0.5 GW tri , 
with 2 years of data taking, and the second one ("LARGE") is a power reactor (10 GW tn , 5 
years data taking) located at a suitable distance from the detector. It is desirable that this 
reactor should not be running during the whole lifetime of the detector, so that the latter can 
also pursue physics goals which require low backgrounds, e.g. the search for Geo-neutrinos, 
supernova relic neutrinos, and proton decay. Of course, no one would build a detector like 
LENA voluntarily in the direct neighborhood of such a power station, since the neutrino 
flux coming from this reactor would dominate all events from other neutrino sources and the 
detector would be effectively unusable for experiments other than long baseline oscillation 
measurements. However, if the reactor is scheduled to be shut down after the first years of 
data taking with the Large Liquid Scintillator Detector or if it is just planned to be built 
but there is enough time left to take data for the other measurements with the detector, 
this will not cause any problems. 

The detector properties and systematical errors are the same as in Section [31 but since 
background sources, in particular Geo-neutrinos, are much more important now, they need 
to be treated in a more sophisticated way. 

In Refs. [16, 46] it has been shown, that Geo-neutrinos have a strong influence in a "large 
reactor - small detector setup", namely SADO. Since we analyze a much larger detector, 
this influence could be stronger here, even if the "product" of the reactor and the detector 
size (i.e. the total exposure) is similar. However, as will be shown, our results are still 
comparable to those obtained by the previous analyses. 

The x 2 function for Section 4, including contributions from Geo-neutrinos in a more accurate 
way, has the following form: 



x 2 



77 [^i( a norm, a det, O reac , Q\J, Oxh, &) — Aj] + 



Aj 
i 



, a norm , a det , a reac , a U , a Th , ^ / 1fl \ 

H j ' 2 1 2 1 2 H T~ + ~~2' 

°"norm ^det a reac °"u a Th °b 

The meaning of the parameters is the same as in Sec. [3J but the background errors <7bk g 
are now split up into the contributions coming from distant nuclear reactors (a reac ), the 
Geo-neutrinos from the uranium decay chain (o"u), and the Geo-neutrinos from the thorium 
decay chain (uxh)- 

We will consider three different situations to show the impact of the Geo-neutrino back- 
ground: 

• No Geo-neutrinos: In this case, Geo-neutrinos are completely absent, and only the 
background from distant nuclear reactors is taken into account. As in Sec. O, it yields 
850 events per year. The uncertainty in the flux normalization is taken to be 2%. 
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• Geo-neutrinos with a 10% uncertainty: Here, Geo-neutrinos are taken into ac- 
count, and the uncertainty in their flux is assumed to be 10%. We use independent 
normalization factors for neutrinos originating from the uranium decay chain and those 
originating from the thorium decay chain (cf. Eq. (I19p ) because the relative abundances 
of these elements in the Earth depend strongly on the geological model [48] . However, 
as central value, we assume the ratio of the thorium and uranium abundances to be 
3.9, as given in [48]. The decay chains of other long-lived radioactive isotopes, e.g. K- 
40, are not relevant for our discussion because they yield anti- neutrino energies below 
the detection threshold of 1.8 MeV. The Geo- neutrino spectra in our simulations are 
taken from [37,38]. The reactor background is the same as in the scenario without 
Geo-neutrinos. 

• Geo-neutrinos with a 100% uncertainty: This scenario is equivalent to the pre- 
vious one, but now the uncertainties in the two Geo-neutrino contributions are taken 
to be 100%. 



4.1 Analytical estimates 



The optimum baseline for a measurement of (9 12 can be estimated analytically if we ne- 
glect backgrounds and systematical errors for simplicity, and perform a total rates analysis, 
neglecting spectral information. Again, we start with the x 2 -function for this scenario: 

2 _ [N (1 - sin 2 26 12 sin 2 A) - N (l - sin 2 29 12 sin 2 A)] 2 
X = iV(l-sin 2 2# 12 sin 2 A) ' (2 ° } 

Here, the barred values are the theoretical predictions, calculated with fit-values for the 
parameters, and the ones without bars come from the observed event rates, calculated with 
the assumed true parameter values from Eq. (j3J). The A's are defined as A = Am 2 1 L/(4i?) 
and A = Am 21 L / (AE) , where L denotes the baseline and E the energy of the incoming 
neutrino. For convenience, we additionally introduce the abbreviations s = sin2#i2 and 
s = sin 2#i2, and write the unoscillated event rate as iV = No/L 2 , where iVo is independent 
of L. Furthermore, for the analytical estimate of the optimal baseline for a measurement 
of the solar mixing angle, we neglect parameter correlations between the solar parameters 
and assume a fixed Am^ = Am^. Now, the x 2 -function becomes: 

^= 717^^ ('-*)'■ <«) 
L z (1 — s z sin A J 

Since we want to measure the amplitude sin 2 26 12 of the oscillations, we can already guess 
that the optimum baseline should correspond to an oscillation phase of approximately tt/2, 
because there, the oscillation has a maximum, and tt/2 is the closest point to the detector 
with this property. With a simple Taylor expansion, one then gets: 

8 mA = s 1 n(f + (A-^))«l-I(A-^) 2 , (22) 

and hence, also up to leading order in the small quantity (A — |J , 

7T \ 2 . 4 , _ „ / . 7T 



sin 2 A » 1 - ( A - - J , sin 4 A«1-2(A--) , (23) 
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which leads to 



L l-« 2 (l-(A-|) ) 



Setting equal to zero, we get: 



2 
L 



70 2 



(A ) -2i(A )s 2 ^ ^ r- = 0. (25) 



£ T 2E y 2 

To lowest order in (A — |), this leads to a baseline of 



1_ 3 2 (l-(A-f) : 



^"^{^f-l^TW)}- (26) 

where i = tan2# 12 . Only the plus-solution is self-consistent with our initial assumption 
A fa tt/2. For E = 4 MeV, it yields a best baseline of approximately 55 km. 

For the determination of Am^, the optimum baseline will be somewhat different, but in 
this case, the calculation is not as simple because L appears both inside and outside the 
sin A terms, so Am^ — Am^ cannot be extracted. Approximations are quite cumbersome 
because of the complex interplay of the different terms, and spectral effects complicate the 
calculation even further. One can however estimate that the best baseline should correspond 
to an oscillation phase A for which small variations of Am^ will have the strongest effect 
on the oscillation probability. It can be read off from Eq. (TJJ) that this will be the case if A 
is an odd integer multiple of 7r/4: A = 7r/4, 37r/4, 57r/4, . . . . Low multiples are favored by 
higher statistics (due to the geometrical factor L~ 2 ), while for higher multiples the effect is 
larger because a variation of Am^ implies a stretching of the sin 2 function in Eq. (jSJ), which 
is proportional to L. This discussion shows that reliable estimates for the optimal baseline 
can only be obtained numerically. 



4.2 Numerical results 



Looking at Figure [31 one can see that the optimum baselines for a measurement of the solar 
mixing angle sin 2 2(9x2 for the SMALL scenario will be around 50 to 70 km if one includes 
Geo-neutrinos, and a bit shorter without Geo-neutrinos (approximately 40-60 km). This is 
in agreement with the optimum baselines obtained for the SADO scenario in Refs. [16, 46] 
(where the technically different situation "large reactor - small detector" is investigated) 
and with our analytical estimates. Geo-neutrinos only perturb the low energy part of the 
spectrum up to around 3.3 MeV, which means that, if they are present, the overall sensitivity 
will be dominated by the high energy part. However, for high energy neutrinos the optimum 
baseline will be larger due to the L/i?-dependence of the oscillation probabilities. For the 
LARGE scenario, all baselines are more or less equivalent, as long as one is far enough away 
from the reactor to be able to resolve the solar oscillation, i.e. the oscillatory disappearance 
"dip" in the spectrum has to be seen. For the determination of Am^, Figure [3] and Tabled 
show that the optimum baseline turns out to be 47.0 km for the SMALL and 53.2 km for the 
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Influence of Geo-v's on #12 (true value: sin 2 2#i2=0.83) 




30 40 50 60 70 80 90 100 
Baseline [km] 



Influence of Geo-v's on Am 2 ol (true value: Am 2 ol =8.2xl0 5 eV 2 ) 




7xl0 -5 



30 40 50 60 70 80 90 100 
Baseline [km] 

Figure 3: Achievable precision for the parameters sin 2 2#i2 (upper panel) and Am^ (lower panel) at the 
90% confidence level as a function of the baseline. The outer curves correspond to the SMALL scenario, 
while the inner ones show the achievable precision withe the LARGE scenario. The true values we have 
chosen are indicated by the horizontal lines. 
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sin^ 26»i2 




Lopt [km] 


events/year 


precision (90% CL) 


SMALL 


no- Geo z/'s 


43.5 


1866 (1026) 


+4.10%/-3.61% 




10% error 


57.5 


2853 (514) 


+5.30%/-5.30% 




100% error 


59.4 


2828 (488) 


+5.42%/-5.42% 


LARGE 


no- Geo z/'s 


40.0 


27572 (26732) 


+0.60%/-0.36% 




10% error 


52.9 


14405 (12065) 


+0.60%/-0.60% 




100% error 


55.0 


13466 (11127) 


+0.72%/-0.60% 




Amy 




Lopt [km] 


events/year 


precision (90% CL) 


SMALL 


no Geo- z/'s 


47.0 


1655 (815) 


+3.05%/-3.54% 




10% error 


76.2 


2729 (389) 


+5.37%/-5.73% 




100% error 


79.7 


2717 (378) 


+5.49%/-5.85% 


LARGE 


no Geo- z/'s 


53.2 


12757 (11917) 


+0.37%/-0.37% 




10% error 


65.0 


11116 (8776) 


+0.37%/-0.37% 




100% error 


70.0 


10589 (8250) 


+0.49%/-0.49% 



Table 2: The achievable precision on the solar parameters sin 2 2#i2 and Am|i at the 90% confidence 
level. Also shown are the optimal baselines and total (signal) event rates per year for the measurements 
in the LARGE and SMALL scenarios. Note that for the LARGE scenario, the exact baseline is not so 
important since one always has extremely high rates. 

LARGE scenario, as long as the Geo-neutrino background is neglected. With the inclusion 
of this background, the best baseline for the SMALL scenario is shifted to 75-80 km, while 
for the LARGE scenario, the impact of the background is less dramatic. 

To get a better overview and to get a clue of how strongly the different scenarios are 
influenced by Geo-neutrinos, we summarize the numerical results in Table [2j Taking into 
account, that our values give the 90% confidence level ranges for the solar parameters, 
one can easily see that already the SMALL scenario could reach values comparable to 
SADO [16,46]. It is important that, for a shorter baseline, the impact of Geo-neutrinos 
broadens the 90% confidence level range quite strongly, even if the flux error is only taken 
to be 10%. Hence, for performing a reactor experiment at such baselines, it is necessary to 
have a good knowledge of the Geo-neutrino fluxes. This goal can be reached for example by 
measuring the Geo-neutrino background before turning on or after turning off the mobile 
reactor. However, for the LARGE scenario, one can clearly see that an amazing precision 
in the measurement of the solar parameters is possible, even without precise knowledge on 
the Geo-neutrino background. One could of course assume, that such a precision could 
be spoiled by the fact that the energy of a single neutrino cannot be determined better 
than about 1 or 2%. However, note that this precision is just the limit for each single 
bin. Since we consider a total of 67 bins the overall precision can indeed be better. More 
precisely, 2% is the absolute uncertainty for each bin, but since information from different 
bins is compared, the relative uncertainty between the bins is the main limiting factor and 
this can, due to the spectral information, be much less than the former one. Hence, the 
precision of our results is realistic for the considered scenarios. 



14 



5 Conclusions 



We have performed analytical and numerical calculations to estimate the potential of reactor 
anti-neutrino disappearance measurements with a Large Liquid Scintillator Detector like 
LENA. For the measurement of the small mixing angle sin 2 26 ) i 3 , we have assumed the 
reactor to be mobile. This allows subsequent measurements in a near and a far position, 
so that many systematical uncertainties are canceled. We have distinguished between two 
different scenarios, a nuclear-powered ship (e.g. an icebreaker or a submarine) and a land- 
based mobile reactor, where the cancellation of systematical uncertainties is different. In 
the case of a nuclear-powered ship, detector related and source related systematical errors 
are eliminated, so that the sin 2 26>i 3 sensitivity follows the statistical limit with increasing 
exposure. In the case of the land-based mobile reactor, the source related flux normalization 
error remains because the reactor needs to be shut down for the displacement. If an exposure 
of 13.5 GW kt yrs is assumed, the nuclear-powered ship scenario can provide a limit of 
sin 2 29 13 < 0.004 at the 90% confidence level, whereas the land-based scenario can achieve 
a limit of sin 2 2# 13 < 0.02 with the same exposure. Furthermore, we have shown that for 
the nuclear-powered ship scenario the two baselines, near and far, should be chosen around 
the 0th oscillation minimum (L « km) and 1st oscillation maximum (L 1.3 km), 
while the land-based scenario yields optimal results with near and far baselines around the 
1st oscillation minimum and 1st oscillation maximum. This is advantageous because the 
detector will be located underground, and the mobile reactor cannot be positioned in the 
direct neighborhood of the detector. 

We have also analyzed the potential of a Large Liquid Scintillator Detector to perform 
precision measurements of the solar parameters sin 2 29^ and Ara^. For this, we have 
compared two scenarios: a small 0.5 GW reactor (SMALL) and a large 10 GW power station 
(LARGE). As has previously been shown in [16], the measurement of the solar parameters 
is strongly influenced by the Geo-neutrino background, therefore we have implemented an 
accurate treatment of this background, and compared its impact for different assumptions 
on the systematical uncertainties. We have shown, that the SMALL scenario would favor 
baselines of approximately 40 to 60 km for the measurement of sin 2 2#i2 and 47 km for the 
measurement of Ara^. The exact choice of the baseline in the LARGE scenario turned out 
to be less crucial, since the excellent statistics provide precise information on the energy 
spectrum. Even with the most conservative assumptions on the systematical uncertainties 
of the Geo-neutrino background, the SMALL scenario can already achieve precisions of 
approximately 5.4% on sin 2 2#i2 and 5.9% on Am^, while the LARGE scenario achieves 
precisions of approximately 0.7% on sin 2 2# 12 and 0.5% on Ara^. This very good accuracy 
is required e.g. to test quark-lepton complementarity at a very high precision level. 
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